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Abstract. - We study a one-dimensional Bose system with repulsive 5-function interaction 
in the presence of an SU(2) intrinsic degree of freedom on the basis of the coordinate Bethe- 
ansatz. The ground state and the low-lying excitations are determined by both numerical and 
analytical methods. It is shown that the ground state is an isospin-ferromagnetic state, and 
the excitations are composed of three elementary particles: holons, antiholons, and isospinons. 
The isospinon is a triplet coupled to the "ferromagnetic" background anti-parallelly. 



Exactly solvable models @,|,|,|,|,|,0,|,|,0,|T^,|T|,|T|,[T4[|T|,[6) play an important role 
in physics, specifically in the investigation of one-dimensional (ID) interacting many-particle 
systems. They have served as a source of inspiration for the understanding of non-perturbative 
phenomena in correlated electronic systems; e. g., the spinon was explicitly characterised on 
the basis of the exact solution of the Hubbard model ||. Among these models, an earlier 
prototypical one is the model of ID bosons with repulsive (5-function interaction, which was 
solved @ by means of the Bethe-ansatz. This method was also applied to solve the problem of 
spin- 1/2 fermions @,[| with S- function interaction; in fact, Yang |9| already suggested a general 
strategy for multi-component systems. Various extensions include the study of electrons on 
a crystalline lattice 0, the generalization to higher symmetries and applications 

to different boundary conditions Jl3|, |l|, [ll| [l6| . Nevertheless, SU(2) bosons with S- function 
interaction, as far as we are aware, have not been studied until now - in contrast, obviuously, to 
the "spin 1/2" fermionic case, i. e. a two-component model with anti-symmetric permutation. 
Recently, however, a two-component Bose gas was created in magnetically trapped 87 Rb by 
rotating the two hyperfine states into each other with the help of a slightly detuned Rabi 
oscillation field [[t7|, [l8| ; and it was noticed |l^] that the ground state of a Bose system can 
be surprisingly different from the scalar Bose system, once the particles acquire an intrinsic 
degree of freedom. 

Therefore, in order to obtain non-perturbative insight into the features of one-dimensional 
SU(2) bosons with repulsive (5-function interaction, we study a model which is integrable. 
Pointing out the connection with the coupled Gross-Pitaevski equation |^,|2lJ], we solve the 
Bethe-ansatz equations for bosons with an SU(2) intrinsic degree of freedom. The ground state 

© EDP Sciences 



2 



properties and the low-lying excitations are studied by both, a numerical calculation and in 
the thermodynamic limit. Unlike a spin-1/2 Fermi system, the ground state of the present 
model is a "ferromagnetic" state, consistent with the result of a variational approach [jTsfl . 
The charge-isospin phase separation is confirmed, and the isospinon is a triplet instead of a 
doublet (the latter being well known for the spinon in the spin-1/2 Fermi system). 

The two-component Bose gas is known to satisfy the coupled Gross-Pitaevski equations: 

dt\i> 2 ) \P G 2 J\i> 2 

where Gb = — T~V^ + Vb(r) + ^ a u a b\ipa\ 2 , P = Ml/2, and a, b = 1,2; CI is the Rabi oscillation 
field frequency, and V(r) the trapping potential. 

We consider the isotropic limit, in which the strengths of interaction between inter-species 
and intra-species are the same {u a b = c); then the model is integrable. Considering further 
the system trapped in a ID ring of length L, and introducing cf>i = (tpi + tp2 ) / V2, 4>2 = ("01 — 
02)/v / 2, we obtain from the above equation (for real f2) the following equivalent Hamiltonian: 

H = J dx[d x <t>* a ■ d x <t> a + CCp* a 4>t4>b^a ~ (-l) B fi##„ 

with summation over repeated indices implied; natural units are adopted for simplicity. The 
fields obey bosonic commutation relations, [<f>* a (x), 4>b{y)] = Sne^^ ab ^( x ~ V ~ nL). The 
Rabi oscillation field contributes a Zeemann-type term. To avoid confusion with conventional 
spins, we denote the generators of the isospin SU(2) by 1; correspondingly = 2T Z . 

The intrinsic degree of freedom can be specified by the eigenvalues of T z , or by isospin up and 
down. 

The Bethe-ansatz equations for two-component bosons are given as follows: 

N M 

e^ L = -H^i(kj - h) J] S_ 1/2 (kj - A„), 

1 = 1 v=l 

N M 

l = -I|M A 7^i)II 5 l( A T-U 

1 = 1 v=l 

where ^p{x) — (x + if3c)/(x — if3c). Equation (|l|) determines the value of the quasi-momenta 
{kj} and the isospin rapidities {A„} for a N — 2M + 1 fold multiplet characterised by the total 
isospin X to t — (N — 2M)/2. These equations are similar to the Bethe-ansatz equation of |22| , 
except for a variation in the exponential. 

Equation (|l|) is obtained as follows. Applying the Hamiltonian to the Hilbert space of 
N particles, and considering its first quantized version on the domain R \ {Py} where Pjj := 
{x\xi — Xj — 0} is the hyperplane defined by the ^-function singularity, we see that only the ./V- 
dimensional Laplacian remains in the Schrodinger operator. Thus TV-dimensional plane waves 
are solutions. We sum up all the plane waves with wave vectors which are are just permutations 
of a definite k = Qti, ■ ■ ■ , /sat) according to the Bethe-ansatz strategy. Integrating the 
Schrodinger equation across the hyperplanes, we obtain S(ki—kj) — [ki — kj — icP]/[ki—kj+ic], 
which connects the wave functions defined on the regions separated by the hyperplanes, and 
S := VS (V stands for the spinor representation of the permutation group Sn) relates the 
coefficients of different plane waves in the same region. The bosonic permutation symmetry 
(instead of the antisymmetry) was imposed when solving for the S'-matrix. Analogous to 
the case of spin-1/2 fermions ||, the periodic boundary condition leads to an eigen-equation 
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for the product of the ^-matrices. As the S'-matrix satisfies the Yang-Baxter equation, the 
quantum inverse scattering method p3[ is applicable. After writing out the fundamental 
commutation relations, and evaluating the eigenvalues of the reference state \u>) = | ft • • • 1)i 
one immediately recognises the differences to the case of spin- 1/2 fermions. For example, 

M0\") = TIM - 6 " ic)/(£ - £, + icJIw), and £>(£)M = ]!;(£ " ~ 6 + "OM, m the 

notion of [p3"[p4j . Consequently, we obtain Eq. (|l|). We note that the Bethe-ansatz strategy 
implies the existence of infinitely many constants of motion, Y) , kj — constant, in addition 

to the usual energy E = Y2i=i kf + $l(N — 2M), and momentum P = Y^f hi. 
Taking the logarithm of Eq. (Q) leads to 



N * M 



l=l v=l 
N M 

2ttJ 7 = ^6_ 1/2 (A 7 - h) + ^9i(A 7 -A„), (2) 



where Qp(x) :— — 2 tan _1 (x//3c); both the quantum numbers Ij and J 7 take integer or half- 
integer values, depending on whether N—M is odd or even. For comparison, the Bethe-ansatz 
equation for spin- 1/2 fermions does not only lack the first summation, but also has an opposite 
sign in the second summation in the first line of Eq. (Q) . The momentum is easily obtained 
from Eq. (§), P = £, *, = (£, I* - J !/ )2 7 r/L. 

It is instructive to analyse Eq. (^|) in the strong and weak coupling regimes. For a strong 
interaction, c — * oo, the wave function vanishes for any x% — Xj, and hence the bosons avoid 
each other like fermions which is in agreement with the discussion of quantum degeneracy 
in trapped ID gases |p5| . On the other hand, in the weak coupling limit, c — » 0, using 
Gi(x) — > — 7rsgn(x), and 9 _ 1/2(2;) — > 7rsgn(a;) for x 1, Eq. (||) becomes 



1=1 v=l 

N M 

X S § n ( A 7 - k l) - ^2 S § n ( A 7 - A ») = 2J 1- ( 3 ) 

i=l !/=l 

The subscript of the isospin rapidity A 7 can be chosen in such a way that J 7 is arranged in 
increasing order; then the second equation of Eq. (||) turns into 

N 

^sgn(A 7 -fc;) = 2 J 7 + 27 - M - 1. (4) 
1=1 

Because | J 7 | < (N — M + l)/2 for a given M and M < N/2 due to the restriction given by the 
Young tableau, the minimum value of the right hand side of Eq. (^) is —N + 2. This requires 
that the smallest ki must be smaller than the smallest X v , otherwise the left hand side would 
be —N. Equation (El) also implies 



N 

E 

1=1 



[sgn(A 7+ i - ki) - sgn(A 7 - fej) = 2( J 7+i - J 7 + 1). (5) 
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Thus, for J 7 +i — J-y — m, there must exist exactly m + 1 solutions of ki satisfying A 7 < ki < 
A 7 +i. Furthermore, from the first equation of Eq. (|^), we obtain 

7T M 2lT 

kj+i - k } ;- ^ [ a EP(kj+i - K) - sgn(fcj- X v ) = - Ij - 1). (6) 
i/=i 

Obviously, for Ij+i — Ij = n, there will be fcj+i — kj = 2mr / L if there is a A 7 such that 
kj < A 7 < kj+i, otherwise kj+i — kj = (n — l)2ir/L. Thus an isospin rapidity of value 
A^ always repels the quasi-momenta away from that value. As a result, an existing A M will 
suppress the density of states in fc-space at the point k = A^. The more isospin rapidities 
there are, the higher the energy will be. Thus the ground state of SU(2) interacting bosons is 
no longer a SU(2) singlet, but an isospin "ferromagnetic" state which differs from the Fermi 
case considerably. 

For N particles, the ground state is characterised by a one-row TV-column Young tableau 
[N], of which the quantum numbers are {Ij} := {— (iV— 1)/2, (JV— 1)/2} and {J 7 } = empty. 
For this state Eq. (||) reduces to the case studied in (7) , but the ground state of the present 
model is an (N + l)-fold multiplet with I 2 = N(N + 2)/4. The density of states per length 
for the ground state is plotted in Fig. |l| (left) for various coupling constants. The "particle" - 
hole (or maybe to be called holon-antiholon) excitation is defined by the quantum numbers 
h = ~{N - l)/2 + S ltjl (for 1 < ii < N), Ij = Ij-i + 1 + S jdl (for j = 2, N - 1), and 
\In\ > (N + l)/2. Figure [2] (left) shows the corresponding excitation spectrum. The isospinon- 
holon excitation is characterised by the Young tableau [N — 1, 1], i. e., M = 1. In comparison 
to those of the ground state, the quantum numbers {Ij} change from half-integer to integer 
or vice versa; accordingly, I\ = —N/2 + S\ i j 1 (for 1 < ji < N + 1), Ij = Ij-i + 1 + Sjj 1 (for 
j = 2, ...,N), while J\ = I\ + n so that I\ < J\ < In- This is an (N — l)-fold multiplet 
with X 2 = N(N — 2)/4. The excitation spectrum is plotted in Fig. g. In fact two branches of 
the quasi-particle excitations were recently observed in a two-component condensate p6[ by 
means of techniques from nonlinear optics. The density of states for J\ = 0, ji = 1 is plotted 
in Fig. |l|. In comparison to the ground state where no isospin rapidity exists, a rift emerges 
at the position of the isospin rapidity for small c that is consistent with our previous analysis 
for weak coupling. 

In the thermodynamic limit, the Bethe-ansatz equations lead to the following integral 
equations for the density of roots and the density of holes, respectively, in quasi-momentum 
and isospin rapidity spaces: 

P {k)+p h {k) = -!-+ / dk'pik^K^k-k 1 ) - [ dX'a(X')K 1/2 (k-X'), 
2tt J q J_ b 

<r(X)+a h (X) = / dk'p(k')K 1/2 {X-k') dAV(A')^i(A-A'), (7) 

J~Q J-B 

where K^(x) — tt^ 1 pc/ (p 2 c 2 + x 2 ). The limits of integration, Q and B, are determined 

to be consistent with p(k)dk = N/L, and J B B cr(X)dX — M/L. It is easy to check by 
Fourier transformation that the state with B = 00 and oh — is an isospin singlet, but it 
is not the ground state. The ground state corresponds to a = ph = in Eq. (^), i. e. it 
is an isospin "ferromagnetic" state, in agreement with the result of mean field theory ]l9[| . 
The two-particle case is a pedagogical example: For the two-body Schrodinger equation in 
the center-of-mass frame, the permutation of particle coordinates equals the parity reflection 
of their relative coordinate. The oscillation theorem in quantum mechanics tells that the 
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spatial wave function without nodes, an even parity solution, yields the lowest energy. If it 
possesses an SU(2) intrinsic degree of freedom, the intrinsic wave function must be symmetric 
(anti-symmetric) to keep the total wave function with the lowest energy being symmetric (anti- 
symmetric). Then the ground state of the Bose system (Fermi system) is of "ferromagnetic" 
( "anti-ferromagnetic" ) character. 

The highly degenerate ferromagnetic ground state, obtained in the case of a vanishing 
Zeemann term, will split up into Zeemann sublevels once the external field is applied. The 
ground state hence becomes a polarized state once the Rabi field, which breaks the SU(2) 
symmetry, is turned on. 

In order to evaluate the excitation energies we put p(k) = po(k) + p±(k)/L (po refers to the 
ground state) . In the presence of the isospin degree of freedom, there will be a holon-isospinon 
excitation, in addition to the holon-antiholon excitation. The latter is created by a hole inside 
the quasi Fermi sea k £ [— hp, kp], and an additional k p outside it, i. e. 



The excitation energy consists of two terms: AE = J k 2 pdk + k 2 = £h{k) + e a (k p ), where the 
holon energy £/,, and the antiholon energy e a (k p ) — —eh{k p ), are given by 



Flipping one isospin corresponds to adding one isospin rapidity to the background of the 
ferromagnetic ground state, which inevitably brings about one hole in the /c-sector. The 
excitation energy AE = J k 2 p\dk is obtained from 



In conclusion we found three elementary quasi-particles: holon, antiholon and isospinon. 
From the asymptotic behaviour of those basic modes, for k, k p , and A tending to kp, we 
find that both the holon-antiholon and holon-isospinon excitations are gapless. The related 
dispersions for finite N are plotted in Fig. ||. Different from the spinons in a Fermi system, 
the isospinon here is a triplet which always couples to the "ferromagnetic" background anti- 
parallelly. Although it is always accompanied by a charge excitation (holon) for single or odd 
number of isospinons, the isospinons can be excited in pairs without exciting the U(l) charge 
mode. Because of the coupling between the charge sector and the isospin sector in Eq. (R), 
both cases bring about changes in the quasi-momentum distribution and hence lead to an 
excitation energy. The charge-isospin phase separation predicted by mean field theory (^?J 
is clearly confirmed in the present case due to the structure of Eq. (Q). The holon and 
antiholon are quasi-particles created in momentum space, while the isospinon behaves like a 
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dark soliton [£8| in the isospin sector that tends to decrease the total isospin X eigenvalue by 
one. 

Considering the experiment on a two-component Bose gas whose transverse excitations are 
frozen out, such that the dynamics becomes essentially one-dimensional pfij ], it is expected 
that the above results can be confirmed through a careful measurement of the excitation 
spectra, which in particular should show the isospin excitations. However, the question has to 
remain open of how much can be learned from our model - which is an integrable one - about 
the phenomenon of Bosc-Einstcin condensation, as e. g. the inclusion of a realistic trapping 
potential into the model should be important. 
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Fig. 1 - The density of state per length in fc-space for the ground state (left) and for the state in 
the presence of one isospin rapidity A by choosing Ji = (right). The distribution gradually changes 
from a "histogram" to a narrow peak for strong to weak coupling, c = 10, 1, 0.1, 0.01. The calculation 
is done for L = 40, N = 40. The left panel is similar to Fig. 2 of [7]. 




Fig. 2 - The holon-antiholon excitation spectrum (left) and holon-isospinon excitation spectrum 
(right), calculated for L = 20, N = 40, and c = 10. 




Fig. 3 - Dispersion relations of the isospinon for different coupling constants (left); the curves from 
top to bottom correspond to c = 1, 10, 20, 40, and 80, respectively. Right: dispersion of antiholon, 
holon, and isospinon, for c = 1 (L = 20, N = 40). 



